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ABSTRACT 

An analysis of the X-ray variability of the low luminosity Seyfert nucleus NGC 4395, based 
on a long XMM-Newton observation, is presented. The power spectrum shows a clear break 
from a flat spectrum (a « 1) to a steeper spectrum (a w 2) at a frequency /t, r = 0.5 — 3.0 x 
10~ 3 Hz, comparable to the highest characteristic frequency found previously in a Seyfert 
galaxy. This extends the measured Mbh — /br values to lower Mbh than previous studies of 
Seyfert galaxies, and is consistent with an inverse scaling of variability frequency with black 
hole mass. The variations observed are among the most violent seen in an AGN to date, with 
the fractional rms amplitude (F val ) exceeding 100 per cent in the softest band. The amplitude 
of the variations seems intrinsically higher in NGC 4395 than most other Seyferts, even after 
accounting for the differences in characteristic frequencies. The origin of this difference is not 
clear, but it is unlikely to be a high accretion rate (L/LEdd <> 20 per cent for NGC 4395). 
The variations clearly follow the linear rms-flux relation, further supporting the idea that this 
is a ubiquitous characteristics of accreting black holes. The variations are highly coherent 
between different energy bands with any frequency-dependent time delay limited to <J 1 per 
cent. 
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1 INTRODUCTION 

The nucleus of the dwarf galaxy NGC 4395 (Filippenko & Sargent 
1989) lies at the edge of the parameters space formed by Active 
Galactic Nuclei (AGN). It is among the nearest (D ~ 4.0 Mpc; 
Thim et al. 2004) and lowest luminosity (LboI ~ 5xl0 40 erg s _1 ; 
Moran et al. 1999; Lira et al. 1999) Seyfert galaxies known, and 
harbours a central black hole with a mass probably in the range 
Mbh ~ 10 4 - 10 5 M (Filippenko & Ho 2003), which compares 
to more luminous Seyfert galaxies with black hole masses typically 
in the range Mbh ~ 10 6 - 10 s M Q (Wandel, Peterson & Malkan 
1999). In addition, it shows some of the strongest X-ray variabil- 
ity of any radio-quiet AGN (Iwasawa et al. 2000; Shih, Iwasawa 
& Fabian 2003) with variations of nearly an order of magnitude in 
a few thousand seconds. Yet in other respects, such as optical and 
ultraviolet emission (and absorption) line properties, it resembles 
its more luminous counterparts (Filippenko & Sargent 1989; Cren- 
shaw et al. 2004). 

The X-ray variations of Seyfert galaxies are characterised by a 
'red noise' power spectrum (Lawrence et al. 1987; M c Hardy 1989), 
approximated by a steep power law at high frequencies [P(f) oc 
f~ a ] with an index a ~ 2 breaking to a flatter index ~ 1 below 
some break frequency /b r (e.g. Edelson & Nandra 1999; Nandra 
& Papadakis 2001; Uttley, M c Hardy & Papadakis 2002; Vaughan, 



Fabian & Nandra 2003; Markowitz et al. 2003; M c Hardy et al. 
2004). The prevailing hypothesis is that /b r scales inversely with 
the mass of the black hole (or, equivalently, the break timescale 
scales with Mbh)- Indeed, the limited data to date are consistent 
with a simple scaling relation for all black holes, from Galactic 
Black Holes (GBHs) in X-ray binaries (Mbh ~ 3-20 M Q ; 
M c Clintock & Remillard 2004) to supermassive black holes in 
AGN (M BH ~ 10 6 - 10 8 M ). If this relation holds then the break 
timescales in Seyferts may be used to infer the mass of the central 
black hole by comparison with the power spectra of well studied 
GBHs such as Cygnus X-l (Belloni & Hasinger 1990; Nowak et 
al. 1999) which contains a Mbh ~ 10 M© black hole (Herrero 
et al. 1995). Uncovering the dependence of the characteristic vari- 
ability frequencies with the fundamental source parameters (mass, 
luminosity, accretion rate, etc.) may help clarify the origin of the 
variability itself. 

However, this simple hypothesis is difficult to test (see dis- 
cussions in Markowitz et al. 2003; M c Hardy et al. 2004 and Pa- 
padakis 2004). GBHs typically possess masses Mbh ~ 10 Mq 
whereas the majority of the well studied Seyfert galaxies have 
masses Af B H > 10 6 M©. Thus the 10 2 - 10 6 M Q region of pa- 
rameter space remains unexplored. NGC 4395, with its unusually 
low mass black hole, falls in this critical region. Previous studies 
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of its X-ray variability (Shih et al. 2003) using ASCA suggested the 
expected power spectral break was present (at /br ~ 3 x 10 -4 Hz), 
but this was poorly determined due to the interrupted sampling of 
the data (caused by the periodic Earth occulations and SAA pas- 
sages of the low Earth orbit). XMM-Newton, with its ~ 2 day orbit 
and high throughput, does not suffer from this problem and has 
been providing the best constraints on the high frequency power 
spectra of Seyferts (Vaughan et al. 2003a; Vaughan & Fabian 2003; 
M c Hardy et al. 2004). The present paper describes a time series 
analysis of the data from a long ~ 100-ks XMM-Newton observa- 
tion of NGC 4395 designed to confirm the break timescale and test 
the Mbh — /br scaling hypothesis. A future paper (Iwasawa et al. 
in prep.) will discuss other aspects of this observation. 



2 XMM-Newton OBSERVATION 

NGC 4395 was observed by XMM-Newton (Jansen et al. 2001) 
over 2003-November-30 and 2003-December-01 for a duration of 
113, 384 s. This paper will concentrate on the results from the EPIC 
pn camera which provided the highest signal-to-noise X-ray data. 

Extraction of science products from the Observation Data 
Files (ODFs) followed standard procedures using the XMM- 
Newton Science Analysis System v5.6.0 (SAS). The data were 
taken with the EPIC cameras in full-frame mode. The field of 
NGC 4395 contains several other bright X-ray sources, although 
the nucleus itself is not confused by these in the EPIC pn image. 
Source data were therefore extracted from the pn image using a 
circular region (of radius 45 arcsec) around the centroid of the nu- 
clear source. Only events corresponding to patterns — 4 (single 
and double pixel events) were used for the pn analysis. Background 
events were extracted from a region on the same CCD unaffected 
by source photons. These showed the background to be relatively 
low and stable throughout the first 90 ks of the observation. During 
the final ~ 15 ks (as the spacecraft approached the radiation belts 
at perigee) the background rate increased dramatically; these data 
were ignored for the present analysis. The final extracted dataset 
contained « 9.85 x 10 4 source counts. The source count rate was 
low enough that pile-up effects were negligible. 

Light curves were extracted from the EPIC pn data in four 
different energy bands: 0.2 - 10.0 keV (full band), 0.2 - 0.7 keV 
(soft band), 0.7 - 2.0 keV (medium band) and 2.0 - 10.0 keV (hard 
band). These were corrected for telemetry drop outs (less than 1 per 
cent of the total time), background subtracted and binned to 50 s 
time resolution. The errors on the light curves were calculated by 
propagating the Poisson noise. The light curves were not corrected 
for the ~ 71 per cent 'live time' of the pn camera (which is only a 
scaling factor). The full band light curve is shown in Fig.Q During 
the XMM-Newton observation the average 2—10 keV flux was 
ss 5.6 x 10 -12 ergs -1 , similar to the flux during the long ASCA 
observation taken in 2001 (« 4.7 x 10~ 12 ergs -1 ; Shih et al. 
2003), indicating the source was in a fairly typical flux state. 



3 TIME SERIES ANALYSIS 

The average source (background subtracted) and background count 
rates are shown in Table Q along with the fractional excess rms 
variability amplitude of the source (F val ; Vaughan et al. 2003b) in 
each energy band. The variability amplitude is remarkably high, 
as can clearly be seen in the light curve (Fig.0. A fractional rms 
in excess of unity is extremely rare and only occurs in the most 




Time (s) 

Figure 1. EPIC pn light curve of NGC 4395 (0.2 - 10 keV with 100 s 
bins). The upper panel shows the light curve with the ordinate plotted on a 
linear scale. The lower panel shows the same time series using a logarithmic 
scale for the ordinate. 



Table 1. Basic properties of the light curves. Mean source and background 
count rates and fractional excess rms variability amplitudes (i*Var) in each 
of the four energy bands (full = 0.2 - 10.0 keV; soft = 0.2 - 0.7 keV; 
medium = 0.7 - 2.0 keV; hard = 2.0 - 10.0 keV). 



Band 


source (ct s 1 ) 


background (ct s 1 ) 


Fvar (%) 


full 


1.11 


1.9 X 10~ 2 


70.4 ± 1.2 


soft 


0.28 


5.9 x 10~ 3 


105.9 ± 1.9 


medium 


0.34 


5.0 x 10~ 3 


83.8 ±1.5 


hard 


0.49 


7.8 x 10~ 3 


43.2 ±0.9 



variable AGN known (NGC 4051, Green, M c Hardy & Done 1999; 
IRAS 13224-3809, Boiler et al. 1997). Thus NGC 4395 clearly 
deserves to be compared to the most variable AGN. The maximum- 
to-minimum variation through the length of the observation is a 
factor ~ 27. 

3.1 Rms-flux correlation 

As can clearly be seen in the light curve, the distribution of fluxes 
seems asymmetric, with a skewed tail of very high fluxes, inter- 
spersed between periods of low and stable flux. (As a crude test, 
examine the light curve upside down. It appears qualitatively dif- 
ferent when viewed this way indicating an asymmetry.) One might 
be tempted to naively interpret the high points as individual events 
('flares' or 'shots') occurring between relatively quiescent periods. 
However, there is a more straightforward interpretation in terms of 
the 'rms-flux' relation known to exist in both Galactic X-ray bi- 
naries (Uttley & M c Hardy 2001) and Seyfert galaxies (Vaughan 
et al. 2003a; Vaughan et al. 2003b; IvTHardy et al. 2004). FigH 
demonstrates this relation operates in NGC 4395. As discussed by 
Uttley, IvTHardy & Vaughan (2004) this is a natural artifact of 
a non-linear transformation of some underlying linear variability 
process. At high fluxes the rms amplitude is proportionally higher, 
hence the light curve appears more erratic, while at low fluxes the 
rms is lower and the light curve smoother. In effect the non-linear 
transformation is a 'stretching' at high fluxes and a 'compression' 
at low fluxes. The effect on the light curve is more apparent in 
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Figure 2. Rms-flux correlation. The flux and noise-subtracted rms were 
calculated every 500-s from the (50-s binned) full band light curve. These 
were then binned as a function of flux such that each bin contains 20 
estimates. Clearly the mean rms is correlated with the mean flux. The dotted 
line shows the best-fitting linear function (\ 2 = 8.09/7 dof) 
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Figure 3. Log-binned periodogram. Upper panel: The histogram shows the 
binned data (TV 20 per bin) with associated error bars. The dotted lines 
shows the expected Poisson noise level. The solid curve marks the best- 
fitting broken power law model (see text). Lower panel: residuals from the 
model fit. 

NGC 4395 simply because of its exceptionally large variability am- 
plitude, which increases the overall effect. In order to illustrate this 
point the lower panel of Fig. Q shows the same light curve with 
a logarithmic flux scale, which effectively transforms the distribu- 
tion of data to a more linear state and eliminates the linear rms-flux 
correlation. (On the logarithmic scale the plot does not look so dif- 
ferent from its upside down counterpart.) Thus the appearance of 
'flaring' and 'quiescent' periods does not require any non-linear 
dynamics but is simply a consequence of the well-established rms- 
flux relation (which is a static non-linear effect). 



3.2 Power spectral analysis 

The power spectrum was examined by calculating the logarithmi- 
cally binned periodogram (following the prescription given in Pa- 
padakis & Lawrence 1993). The result is shown in Fig. [3] which 



clearly reveals the red noise spectrum spanning at least two decades 
of temporal frequency. 

In order to compare the spectrum with those of other Seyferts, 
as well as GBHs, and to determine the interesting power spectral 
parameters, simple models were fitted to the binned periodogram. 
The Poisson noise in the light curve added a constant 'background' 
level to the spectrum, which was accounted for in the fitting by 
including a constant at the appropriate level in the model (see Ap- 
pendix A of Vaughan et al. 2003b). It is important to note that this 
observation of NGC 4395 is free from spectral distortions such as 
'aliasing' and 'red noise leak' that commonly affect Seyfert obser- 
vations. This was demonstrated by comparison with Monte Carlo 
simulations (see below), and implies that direct fitting of an analyt- 
ical model to the binned data is a perfectly reasonable procedure. 

As a first test a simple power law was fitted to the binned peri- 
odogram. This gave a rather poor fit (x 2 = 91 .3/43 dof for the full 
band data). Allowing for a break in the power law gave a substantial 
improvement (Ax 2 = 14.9 for 2 additional free parameters), sig- 
nificant at > 97 per cent in an F-test. The best-fitting broken power 
law model and the data-model residuals are shown in Fig. [3] The 
free parameters were the power law index below the break («i), 
the index above the break (02), the break frequency (/br) and the 
normalisation. Table [2] summarises the results of fitting this model 
to the data from all four energy bands. The low frequency slope is 
consistent with that obtained from an analysis of the long ASCA ob- 
servation by Shih et al. (2003; ax = 0.98±0. 18), but the break fre- 
quency is marginally inconsistent with this previous analysis. This 
may be partly due to the worse sampling (resulting in some spec- 
tral distortion) for the ASCA data, but could also indicate the break 
frequency genuinely changed between the two observations. 

The broken power law model provided a formally unaccept- 
able fit to the full band data, but the residuals (Fig. [3} did not show 
any obvious cause for alarm. The soft and medium band data were 
reasonably well fitted using this model, while the hard band re- 
mained poorly fitted. Part of this may be due to the model having 
the wrong shape break (e.g. too sharp), although fits with smoothly 
bending power laws (see M c Hardy et al. 2004) gave very similar 
quality fits (x 2 = 76.39/41 dof for the full band data). Regard- 
less of whether the shape of the break is more complicated than 
used in these simple models, the existence of a break is clearly 
demonstrated. In all the bands the best-fitting model breaks from 
at\ ~ 1.2 to 02 ~ 2.0 with the break frequency constrained to 
lie in the range / br « 0.5 - 3.0 x 10~ 3 Hz. Fig. H shows the 
change in Xmin as eacn parameter is varied. Using the smoothly 
bending power law model predicted a slightly lower break fre- 
quency (/br = 0.81q g x 10 -3 Hz), compared to the sharply 
breaking model, together with a slightly flatter low frequency slope 
ai w 0.7. 

The Monte Carlo fitting procedure discussed by Vaughan et 
al. (2003a; see also Green et al. 1999 and Uttley et al. 2002) was 
also used to fit these data. This takes into account any power spec- 
tral distortion due to the time sampling of the data. However, in the 
case of NGC 4395 there were found to be no differences in best- 
fitting model parameters between the simple (analytical) and robust 
(Monte Carlo) fitting techniques. There are two reasons for this. 
The first is that aliasing is negligible due to the fact the data are con- 
tiguously binned from an (effectively) continuous photon counting 
signal (van der Klis 1989). The second is that red noise leak was not 
present because the power spectrum was flat ot\ ~ 1.2 at the low- 
est measured frequencies; leakage is only substantial when the low 
frequency power spectrum remains steep (ai > 1.5; Uttley et al. 
2002). Given that the binned periodogram is free from distortion, 
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Figure 4. Parameter uncertainties. Variation in \ 2 f° r the broken power 
law model as the three parameters (111,02 and /br) ale varied. The dotted 
lines show the A% 2 = 2.71, 6.63 bounds, which delineate the 90 and 99 
per cent confidence bounds on the individual parameters. 



Table 2. Results of fitting the periodograms derived from the four en- 
ergy bands (full = 0.2 - 10.0 keV; soft = 0.2 - 0.7 keV; medium 
= 0.7 - 2.0 keV; hard = 2.0 - 10.0 keV) with the broken power-law 
model (breaking from a slope of ct\ below /b r to «2 above). Errors on 
the model parameters correspond to a 90 per cent confidence level for one 
interesting parameter (i.e. Ax 2 = 2.71). 
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the simple fitting was used over the Monte Carlo method because it 
is substantially faster and allowed for a much larger grid of model 
parameters to be computed in a reasonable amount of processing 
time. (For most other Seyferts, where the power spectrum remains 
steep at the lowest frequencies probed by XMM-Newton, the Monte 
Carlo method should be used.) 
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Figure 5. Coherence between energy bands. The three panels show, for 
each combination of the three energy bands, the coherence between the 
variations as a function of frequency. The coherence was calculated using 
the prescription of Vaughan & Nowak (1997) and binned into logarithmic 
frequency intervals ensuring that each bin contained 10 cross-spectrum 
estimates. Clearly at low frequencies the coherence is high > 0.9 in be- 
tween all the bands. At the highest frequencies the affect of Poisson noise 
makes it difficult to measure the intrinsic coherence. 



Fig. [5] shows the coherence between the three different en- 
ergy bands for NGC 4395. Clearly at low frequencies the variations 
are highly coherent (~ 1.0), while at high frequencies the coher- 
ence becomes increasingly difficult to measure due to the relatively 
stronger contribution from Poisson noise. The coherence was cal- 
culated using the method of Vaughan & Nowak (1997) which cor- 
rects for the effect of Poisson noise. There is a slight tendency for 
the coherence to fall off at higher frequencies in the soft-to-hard 
comparison (lower panel of the figure) which would be consistent 
with that seen in other Seyfert 1 galaxies (e.g. Papadakis et al. 2001; 
Vaughan et al. 2003a; M c Hardy et al. 2004). Fig.[6]shows the cor- 
responding time lag spectra. It is clear that there are no signifi- 
cant time delays between the bands, even at the lower frequencies, 
with an upper limit of ~ 100 s at / ~ 10~ 4 Hz (T ~ 10 4 s), 
which corresponds to a < 1 per cent delay. Unsurprisingly the 
cross-correlation functions (computed using the Discrete Correla- 
tion Function of Edelson & Krolik 1988) were reasonably symmet- 
ric and peaked at zero lag. 



3.3 Cross-spectral analysis 

The cross-spectrum is closely related to the power spectrum. It of- 
fers a comparison between two simultaneously sampled time series 
(such as X-ray light curves in different energy bands) in frequency 
space. The modulus of the cross-spectrum is used to define the co- 
herence (Vaughan & Nowak 1997), a measure of the degree of cor- 
relation between variations in the two time series. If one dataset 
is simply a delayed or smoothed copy of the other, the coherence 
will be unity. If the two datasets reflect totally independent varia- 
tions the coherence will be zero. The cross-spectrum can also be 
used to measure the phase lag, or the time delay between the two 
time series (assuming they are reasonable coherent). The coherence 
and phase lags of GBHs have been discussed by e.g. Nowak et al. 
(1999) and examined in Seyfert galaxies by Papadakis, Nandra & 
Kazanas (2001), Vaughan et al. (2003a), Vaughan & Fabian (2002) 
and M c Hardy et al. (2004). 



4 DISCUSSION 

4.1 Summary of results 

The high frequency variability properties of NGC 4395 have been 
examined using a long (~ 100 ks) XMM-Newton observation. The 
variations clearly followed the linear rms-flux relation previously 
observed in both GBHs and other Seyfert galaxies, further support- 
ing the idea that this is a universal property of the X-ray variations 
from accreting black holes. The high frequency power spectrum 
was steep a ~ 2 above a break frequency /b r ~ 2 x 10 -3 Hz, 
below which it broke to a slope a ~ 1. The break frequency is the 
second highest recorded to date for an AGN (the record is presently 
held by Ark 564; Papadakis et al. 2002). The cross-spectrum be- 
tween different energy bands shows the variations at frequencies 
< /br were highly coherent and any frequency dependent time 
delay is constrained to be shorter than r <; 0.01/ _1 .In other 
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Figure 6. Time delay between energy bands. The three panels show, for 
each combination of the three energy bands, the time delay between the 
variations as a function of frequency. The data were binned into logarith- 
mic frequency intervals ensuring each bin contained 10 cross-spectrum 
estimates. There are no clearly detected delays. 
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Figure 7. The Mbh — /br relation for 11 Seyfert galaxies including 
NGC 4395. The masses are from reverberation mapping experiments ex- 
cept for the four objects marked using dotted error bars. Also shown are 
typical break frequencies for Cyg X-l in both its low/hard and high/soft 
states. The dotted lines show example A/bh °c V/br relations consistent 
with the Cyg X-l points. 



Seyferts, such as NGC 7469 (Papadakis et al. 2001), MCG-6-30-15 
(Vaughan et al. 2003a) and NGC 4051 (M c Hardy et al. 2004), the 
delays are typically ~ 1 per cent. The limits of the time delays for 
NGC 4395 therefore did not rule out typical time delay character- 
istics. Although the variations in these relatively broad-band light 
curves were highly coherent, there is a strong energy-dependence 
to the amplitude of the variations. This will be discussed in a forth- 
coming paper (Iwasawa et al. in prep). 



4.2 The M BH /br relation 

Fig. Q shows the Mbh — /br relation for ten Seyfert galaxies 
using measurements from the literature (see also Markowitz et 
al. 2003; M c Hardy et al. 2004; Papadakis 2004), together with 
the new break frequency determined for NGC 4395. The power 
spectral break frequencies from the other ten Seyferts were col- 
lected from the following papers: Ark 564 (Papadakis et al. 
2002); NGC 4051 (M c Hardy et al. 2004); NGC 3516, NGC 3783, 
NGC 4151, NGC 5548, Fairall 9 (Markowitz et al. 2003); Mrk 766 
(Vaughan & Fabian 2003); MCG-6-30-15 (Vaughan et al. 2003a) 
and NGC 5506 (Uttley et al. 2002). Black hole mass estimates 
were taken from the reverberation mapping studies where possible: 
Ark 564 (upper limit from Collier et al. 2002); NGC 4051 (Shem- 
mer et al. 2003); NGC 3516 (Wandel et al. 1999); NGC 3783, 
NGC 4151, NGC 5548, Fairall 9 (Kaspi et al. 2000). For NGC 5506 
the mass was estimated using the A/bh — cr» relation (section 2.3 
of Papadakis 2004). For Mrk 766 and MCG-6-30-15 the mass es- 
timates are very uncertain. Based on its bulge luminosity, Uttley et 
al. (2002) argued for A/ B h ~ 10 6 in MCG-6-30-15. Due to the 
ambiguities in this method the error bars were taken to span an or- 
der of magnitude in each direction. For Mrk 766, Wandel (2002) 
estimated A/bh ~ 10 7 Mq based on the optical luminosity and 
line widths, while Woo & Urry (2002) estimated ~ 3.5 x 10 6 M 
using the same reasoning. For the purpose of comparison a value of 
6.8 x 10 6 Mq (the mean of these two estimates) was used, again as- 
suming an order of magnitude uncertainty in either direction. Also 
shown are the typical frequencies of the high frequency break in 
Cyg X-l in both the low/hard and high/soft states. 



For NGC 4395 the break frequency / bf = 1.9t?;f x 10~ 3 Hz 
was used together with a mass estimate of Mbh = 10 4 — 10 5 Mq. 
This is the 'best guess' of Filippenko & Ho (2003) consistent with 
the available estimates: they measured a firm upper limit on the 
mass of <C 6.2 x 10 6 M from stellar velocity measurements. 
Based on photoionisation modelling of optical emission lines, 
Kraemer et al. (1999) derived an estimate of ~ 1.2 x 10 5 M Q , 
while Filippenko & Ho (2003) estimated ~ 1.3 x 10 4 M based 
on the relationship between optical luminosity and broad line re- 
gion size. These are consistent with the mass estimated from the 
Mbh - a, relation (Tremaine et al. 2002): ~ 6.6 x 10 4 Mq. It is 
clear that, given the significant uncertainties in the mass estimates, 
there is no glaring inconsistency between the new NGC 4395 mea- 
surement and the hypothesised Mbh oc f^ r relation. However, the 
details (such as whether the index of the relation is —1 and which, 
if any, state of Cyg X-l provides the best comparison) remain to be 
investigated as more Seyfert power spectra are measured. 

Hayashida et al. (1998) estimated the black hole masses for 
Seyferts by directly comparing their X-ray power spectra with that 
of Cyg X-l. They measured the frequency f c at which the power 
spectrum reaches the level f c x P(fc) = 10~ 3 and compared 
this to the corresponding frequency in Cyg X-l. The motivation 
for using this frequency, instead of the break frequency, is that the 
break is often harder to measure and (in certain states at least) is 
known to change frequency in Cyg X-l. There are a number of 
underlying assumptions to this method. The first is that both the 
slope and the (relative) normalisation of the high frequency power 
spectra are similar for Seyferts and Cyg X-l, but the frequencies 
scale inversely with the black hole mass. Implicit in this is the as- 
sumption that we are using comparable energy bands in the two 
types of source (the power spectrum is energy dependent). Fur- 
thermore, the original value of f c for Cyg X-l used by Hayashida 
et al. (1998) was taken from an observation of the source in the 
low/hard state. This method therefore also assumes that Seyferts 
most closely resemble GBHs in the low/hard state, which may not 
be the case (Vaughan et al. 2003a; M c Hardy et al. 2004). In the 
case of NGC 4395 the frequency f c is not directly observed; the 
power spectrum lies well above / x P(f ) = 10 -3 at all measured 
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frequencies. Extrapolating the best-fitting broken power law model 
to higher frequencies gives an estimate of f c ~ 7.2 x 10 -2 Hz 
(after accounting for the factor of two difference in the normalisa- 
tion used by Hayashida et al. 1998). Scaling from the correspond- 
ing frequency in Cyg X- 1 , as defined by equation 4 of Hayashida 
et al. (1998) gives an estimate for M BH ~ 6 x 10 4 M Q . This is 
comparable to the estimates made by the other methods although 
if, as discussed below, the normalisation of the power spectrum of 
NGC 4395 is unusually high, this method will not be reliable. 

4.3 A universal power spectrum? 

The early work on the high frequency power spectra of Seyfert 
galaxies found in favour of a universal power spectrum (Lawrence 
& Papadakis 1993; Green et al. 1993). More recent work (Uttley et 
al. 2002; Markowitz et al. 2003) suggests that while most Seyferts 
exhibit a similarly steep (a « 2) power spectrum at high frequen- 
cies, the slope below the break may not be universal, and some ob- 
jects may show a second break (NGC 3783; Markowitz et al. 2003) 
where others clearly do not (NGC 4051; M c Hardy et al. 2004). 

The extraordinary amplitude of the variations in NGC 4395 
suggest that the normalisation of the power spectrum may also vary 
quite substantially from object to object. A fractional rms in ex- 
cess of Fvar > 1 has been noted in only two other Seyferts prior 
to NGC 4395 (these are NGC 4051, Green et al. 1999 and IRAS 
13224-3809, Boiler et al. 1997). NGC 4051 is also thought to con- 
tain a relatively low mass black hole (Shemmer et al. 2003). At 
first sight it remains possible that the high variability amplitudes 
in both NGC 4051 and NGC 4395 could be due to their power 
spectra being shifted to higher frequencies compared to Seyferts 
with more massive black holes (and, since the spectra are red, more 
variance is shifted into the observable frequency range) without the 
need for an intrinsically higher power spectrum normalisation. The 
power at the break multiplied by the break frequency, fh T Pbr, is 
indicative of the total power in the spectrum. For NGC 4395 this 
is ~ 5.4 x 10~ 2 in units of (rms/mean) 2 . The six Seyferts stud- 
ied by Markowitz et al. (2003) have values in the narrow range 
/brPbr ~ 1.6 - 3.0 x 10" 2 , all lower than NGC 4395. Papadakis 
(2004) estimated a 'universal' value of /br-Pbr ~ 1.7 X 10~ 2 for 
Seyfert galaxies (called PSD amp ). This means that, even after ac- 
counting for the different regions of the power spectrum observed 
in different objects, the normalisation of the power spectrum of 
NGC 4395 is intrinsically high. 

What could increase the overall amplitude of the variability 
(at all frequencies)? A plausible suggestion is that an enhanced 
accretion rate (relative to Eddington) somehow provides a more 
unstable accretion flow and hence stronger fluctuations. The ac- 
cretion rate, relative to the Eddington limit, can be estimated for 
NGC 4395 using the available black hole mass estimates, together 
with the bolometric luminosity. Moran et al. (1999) and Lira et al. 
(1999) estimate similar values of LboI ~ 5 x 10 40 erg s _1 . As- 
suming a mass in the range A/bh = 10 4 — 10 5 Mq this gives 
an Eddington fraction of L/Z/Edd ~ 0.4 — 4 per cent. However, 
the mean unabsorbed 2 — 10 keV luminosity during the XMM- 
Newton observation was L2-10 ~ 1-5 x 10 40 erg s _1 . Assuming 
LboI > lOI/2-io would give I/Boi > 1.5 x 10 41 erg s _1 and this, 
less conservative, luminosity estimate gives an Eddington fraction 
of L/L Edd ~ 1.2 - 12 per cent 1 . Thus NGC 4395 is most likely 

1 This assumes the Bolometric luminosity to be related to the mean X-ray 
luminosity by Lb q i > 101/2—10' If however the peak X-ray luminosity 



accreting at a 'normal' rate; super-Eddington accretion seems un- 
likely to be the origin of the unusually large variability amplitude. 
Other possibilities include the black hole spin parameter, a, and 
the inclination of the putative accretion disc-black hole system, i. 
Unfortunately, these are more difficult to test. 
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